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In this letter we examine the production channels for the scalar or pseudoscalar Higgs plus two 
jets at the CERN Large Hadron Collider (LHC). We identify possible signals for distinguishing 
between a scalar and a pseudoscalar Higgs boson. 
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I. INTRODUCTION 

The Higgs mechanism is responsible for electroweak symmetry breaking in the Standard Model (SM) . The experi- 
mental lower limit on the Higgs mass is approximately 114 GeVjlJ. There are many models that contain more than 
one Higgs boson in various numbers of doublets. In the Minimal Supersymmetric Standard Model (MSSM) there 
are two Higgs doublets that give five physical Higgs bosons: two neutral (H±, H®), two charged H^, and one neutral 
pseudoscalar A (for review see ||). In the MSSM the mass limits change slightly with the lightest of the two neutral 
scalars H® (afterwards reffered to as simply H) having a mass greater that about 91 GeV and the pseudoscalar being 
more massive than roughly 92 GeVQ. 

Finding one or more Higgs bosons is the top priority of high energy physics programs around the world. A subset 
of Higgs bosons in some doublet models may be experimentally difficult to distinguish. The characteristics of the 
scalar H and pseudoscalar A Higgs boson within the MSSM are of particular interest. 

We study the production of both a scalar and pseudoscalar Higgs in association with two jets in hadron collisions. 
At the LHC the primary processes that produce a Higgs plus two jets are gg — > ggH and qg — > qgH , accounting 
for approximately 60%(40%) of the total cross-section respectively. The same is true for the production of the 
pseudoscalar. Other channels that contribute to the total cross-section include qq — ► ggH and qq — > qqH, although 
these channels have been shown to add very little to the total cross-section. In the following calculations, only the 
two dominant channels were considered as the other channels are negligible. 

Total cross-sections of the scalar and pseudoscalar plus two jets exist H, 0, |(J at the lowest order. Total cross-sections 
for the inclusive production have been calculated at NLO for the scalar |7], g and for the pseudoscalar Q and at NNLO 
for the scalar |h], O, [lj, |l3[ Q and for the pseudoscalar (ll| |l6| . If we define the K-factor to be the ratio of the higher 
order cross-section to the lowest order, the rate increase at the LHC at NNLO for the scalar inclusive processes [fl2"| 
was reported to be K NNLO (pp — * H + X) = 2 — 2.2 and for pseudoscalar the K-factor Q can be determined to be 
K NNLO (pp -> A + X) = 2 - 2.3 in the mass range M H ,A = 100 - 200 GeV. The total cross-section and the differential 
cross-section for a scalar Higgs plus one jet has been calculated by || [h], |ll|, [TtJ [l^] and the total rate was also shown 
to increase substantially. The NLO corrections to pseudoscalar plus one jet have not yet been computed. In all of the 
processes cited above the rates increased by comparable amounts. We expect our estimates of the Higgs plus two jets 
rates to be conservative, however, since our proposed observable is normalized to the cross-section, we do not expect 
major changes to occur in our analysis at higher orders. 

In this letter, we propose a technique for distinguishing between a scalar and a pseudoscalar Higgs when produced in 
association with two jets by means of a splitting that occurs in a specific integrated operator moment. This distinction 
is important both experimentally and theoretically in order to separate the two kinds of events and understand the 
properties of these particles which would otherwise be very difficult due to the similarity in their physical observables. 

II. EFFECTIVE LAGRANGIAN 

We work in the limit that the top quark is much heavier that the Higgs bosonji], ||, ||, |l^, |2(], integrating out 
the top quark and neglecting all the other quarks that would normally appear in the loop diagrams. This has been 



' bficld@ic.sunysb.edu 



2 



shown to be an excellent approximation and remains very good even when the Higgs mass is heavier that the mass of 
the top quark. In general this approximation is considered to be a good one when My ^ < 2m t . We consider Higgs 
bosons lighter than 200 GeV. The effective Lagrangian used in the scalar case is defined as 

£f ff = -± 9H HG^G a ^ (I) 

where gn = Q. s /3ttv. is the field-strength tensor for the gluons. The vacuum expectation value (vev) of the Higgs 
field is determined in the usual way as v 2 — (\/2Gf)~ 1 and is numerically equal to approximately 246 GeV. For the 
pseudoscalar case we let the Higgs couple to the quarks with a 75 and the effective Lagrangian|2jcan be written as 

Cf B = \g A AG%G a ^ (2) 

where g A = a s /2itv. Here G" = l/2e /i!/pcr Gp . is the dual of the gluon field-strength tensor. 

This effective Lagrangian generates a scalar Higgs coupling to two, three, and four gluons or a pseudoscalar coupling 
to two or three gluons. The four gluon coupling to a pseudoscalar vertex vanishes via the Jacobi identity as it is 
proportional to a completely antisymmetric combination of structure constants. The Feynman rules for these effective 
theories can be found in Q (for the scalar) and || (for the pseudoscalar). 



III. OBSERVABLES AND MOMENTS 

We present our results for the LHC with \/~S = 14 TeV. We have used the CTEQ6L parton distribution functions p^] 
with Ag° = 226 MeV with a one-loop running of a s for consistency with a value of a s (Mz) = 0.137. The transverse 
momentum (p t ) was constrained to be more than 25 GeV for the Higgs and each of the two jets. Also the rapidity 
was constrained to be \y\ < 2.5 for all the outgoing particles. The separation of the jets was restricted to be 
AR = yf&cj? + A77 2 > 0.7. 

The total cross-section of these two channels are shown in Fig. [l]. These cross-sections agree exactly with those in 
the literature |i| once the problems with the effective coupling constants are remedied. When plotted in this linear 
fashion it is interesting to note the differences in the dependence of the cross-sections on the mass of the Higgs boson. 
Both total cross-sections loose more than two-thirds of their value from 100 — 200 GeV and appear in the approximate 
ratio of {gii/9A) 2 =4/9 due to the similarity in their matrix elements. 

Fig. U shows the normalized transverse momentum spectrum of both the production channels. The pseudoscalar 
Higgs pt spectrum was displaced down by 10% to allow the two curves to be distinguished. If this had not been done, 
the curves would lie virtually on top of one another. Fig. ^ shows the center-of-momentum angle between the Higgs 
and the highest p t jet for the two reactions. This shows what would be expected naively, that the Higgs prefers to 
come out back-to-back with the highest p t jet. Once again, the pseudoscalar curve has been scaled down by 20% to 
allow both curves to be seen clearly. No significant differences between these curves were found. 

The authors of [gij presented a technique for determining the CP nature of the Higgs boson in UH production 
based on certain weighted moments of the cross-section. The cross-section integral was weighted by operators Ocp- 
The six operators presented in |24j are scalar and cross products of the momentum of the outgoing particles (in this 
case the massive top quarks). We propose using the same test for the massless quarks and gluons that make up the 
jets. All of these weighted moments were examined as well as some novel ones and the only operator from these sets 
that produced a significant difference between the scalar and the pseudoscalar signals was the operator |24j 

= (pi x z) ■ (ga x z) 
1 101 x2)-(P2 xz)| 
when it was integrated and normalized as prescribed below 

a[0 CP ] = i / OcpdadPS (4) 

where p\ and P2 are the momentum of the two jets and z is the axis of the beam. The a\ operator is sensitive to the 
cosine of the angle between the transverse momentum vectors of the two jets. Distinguishing between the two jets is 
not important as this moment is invariant under 1 <-> 2. Another combination of momentum in the above equations 
that was considered was to use the moment operators presented in p4| with pi = Pmggs an d P2 the momentum of the 
highest p t jet. However, this yielded no differences in the integrated moments making this definition of little use for 
these channels. 
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FIG. 1: Total cross-sections for scalar and pseudoscalar Higgs plus two jets. These curves are for the LHC with the cuts 
described in the text. 



Fig. |j shows the results of this integration as a function of the Higgs mass. If we consider a conservative estimate 
of 100 fb _1 of integrated luminosity at the LHC and take a branching ratio of approximately 10~ 3 as an order of 
magnitude for the decay of the Higgs to a pair of photons, then the a[ai] observable will have a statistical uncertainty 
of about 5%, making these two signals distinguishable at all mass scales. With this conservative estimate on the 
integrated luminosity we would expect to see about 600 scalar events and 1000 pseudoscalar events for a Higgs mass 
of 120 GeV in this channel for this Higgs decay. These numbers are supported by a more detailed analysis using the 
actual branching ratios calculated using hdecay |2jJ] in the tan/? = 1 limit for the pseudoscalar. 

This integrated moment showed a modest (30%) splitting at all Higgs mass scales from 100 — 200 GeV. The 
pseudoscalar does not show much mass dependence. However, the scalar integrated moment rises slightly with 
increasing Higgs mass. This effect might also be useful as another method for constraining the mass of the scalar 
Higgs boson. The splitting in Fig. |] helps to remove the problems created by the degeneracy in the physical observable 
of the scalar and the pseudoscalar. If the two signals could not be separated, the doublet structure of the model would 
not be easily measured. In the case of the MSSM this would mean that part of the supersymmetric signal might be 
lost or the mass of the scalar Higgs may be determined incorrectly if the pseudoscalar events were wrongly identified 
as scalar events. 

Separating the two signals is theoretically intriguing because it appears to be one of the only ways to predict a 
difference between the scalar and pseudoscalar events of this nature at the LHC by means other than the magnitude 
of their cross-section. This is also interesting experimentally as it leads to the possibility of separating the two kinds 
of Higgs events with the added bonus that the z momentum is not needed in this analysis. 



IV. CONCLUSIONS 



The production channels of the scalar or pseudoscalar Higgs plus two jets were found to have many similarities in 
their physical observables and one important difference in the integrated moment a[ai]. This may help to reduce the 
difficulty in distinguishing between the two types of events at the LHC. The most important aspect of separating the 
two signals is to make sure that the doublet structure (the supersymmetric signal in the case of the MSSM) is not 
lost because of its small cross-section and its similarity to the scalar Higgs with respect to its physical observables or 
wrongly determining the mass of the scalar Higgs by misidentifying pseudoscalar events as scalar events. The proposed 
technique presented in this letter may enable these two signals to be separated after a full detector simulation is 
preformed. 
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FIG. 2: Normalized transverse momentum spectrum of the scalar or pseudoscalar Higgs production channels plus two jets. 
The Higgs mass for both the scalar and pseudoscalar is 120 GeV. Note that the pseudoscalar Higgs has been displaced down 
by 10% to allow the two curved to be distinguished. These curves are for the LHC with the cuts described in the text. 
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FIG. 3: Normalized opening angle in the center-of-momcntum frame between the Higgs and the highest pt jet for the scalar 
and pseudoscalar Higgs production channels plus two jets. The Higgs mass for both the scalar and pseudoscalar is 120 GeV. 
Note that the pseudoscalar Higgs has been displaced down by 20% to allow the two curved to be distinguished. These curves 
are for the LHC with the cuts described in the text. 
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FIG. 4: Normalized integrated moment a[ai] for the scalar and pseudoscalar Higgs plus two jets. The center curve shows what 
the observable would look like if both of the processes were to be measured at the same time with degenerate masses. The 
splitting between the two different production channels is clear at all mass scales with a statistical uncertainty of about 5%. 
These curves are for the LHC with the cuts described in the text. 



APPENDIX A: DIFFERENCES IN THE AMPLITUDES 



It turns out that the differences in the scalar (or pseudoscalar) plus two jets amplitudes squared were very small. 
The differences will be presented using the helicity basis presented in ||, || to make for the most compact matrix 
elements squared. These matrix elements have been found to be in exact analytic agreement with the four dimensional 
matrix elements presented in jl7|]. We identify the momentum as follows (where X should be considered the Higgs 
for the process in question, playing the part of either the scalar or the pseudoscalar). All the momenta are outgoing. 



lipi) + q(p2) -> g(-Pa) + g(-Pi) + x(-p 5 ) 
sOi) + gfo) -> g(-pa) + g(-Pi) + x(-p 5 ) 
q{pi) + q(p2) -> q{-pa) + q{-Pi) + x(-p 5 ). 



(Al) 
(A2) 
(A3) 



In the following we define 5 a b = (p a +Pt>) 2 = 2p a ■ px>- Color factors have been included in the expression for the 
qqggH(A) and qqqqH(A) channels as they affect the terms differently but not in the expression for the ggggH(A) 
channel as there is one overall color factor for all the matrix elements squared. Here N is the number of colors. Color 
and spin averages have not been included nor have any coupling constants. 

For the qqggH(A) channel the difference in the scalar minus the pseudoscalar amplitude squared was 15 terms 
out of 626. Setting the color factors to match those presented in 0, Co = (N 2 — 1)/N and C K = (N 2 - 1)N the 
difference was found to be 
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For the ggggH(A) channel the difference in the scalar minus the pseudoscalar was 16 terms out of 2761. The overall 
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color factor is N 2 (N 2 — 1). The difference was 
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Finally, there are two cases for the qqqqH(A) amplitude squared. If there are identical quarks allowed in the 
scattering process (qqqqH (A)) then there are two diagrams that contribute. The color factors here are Ca = N and 
Cf — (N 2 — l)/2N. The difference in the amplitudes squared is 19 out of 39 terms and is equal to 



1^1™* - \M\ 2 q ^ qqA = AC A C F [2 
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If a different quark pair is created (qqq' q' H (A)) , then the difference is smaller as only one diagram is needed for 
the amplitude. Here 6 out of 10 terms survive and are equal to 
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It should also be noted that all these differences are invariant under 1 <-> 2. 
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